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Liischer's finite size mass shift formula in a periodic finite volume, involving forward 
scattering amplitudes in the infinite volume, is revisited for the two stable distin- 
guishable particle system. The generalized mass shift formulae for the boson and 
fermion are derived in the boson-boson and fermion-boson systems, respectively. 
-^ \ The nucleon mass shift formula is given in the nucleon-pion system and the relation 

to the computation within chiral perturbation theory is discussed. 
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1 Introduction 

Nowadays the control of finite volume effects in lattice QCD simulations with 
dynamical fermions becomes a more important issue in order to determine 
the hadron spectrum precisely. Applications of chiral perturbation theory 
(ChPT) [1] to the measured spectrum on the lattice have been done not only 
to achieve the chiral extrapolation [2,3], but also to find its finite volume 
dependence towards the thermodynamic limit [4,5,6,7,8]. 

In this context, Liischer's formula, relating the mass shift in finite volume 
with periodic boundary conditions to forward elastic scattering amplitudes in 
infinite volume, provides us with an elegant tool for such a purpose [9]. In 
the latter reference Liischer presented a rigorous proof of his formula for the 
case of a self- interacting boson to all orders in perturbation theory [10]. By 
writing the difference of the masses between in finite and infinite volumes as 
Am(L) = M(L) — m with the spatial box of the size L 3 , the formula is given 
by 
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where F(u) is the forward scattering amplitude [y = iq Q denotes the crossing 
variable), which has poles at v — ±z/g with v B = m/2. The coupling A is then 
defined by the relation 
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One may ask if one can immediately generalize this formula to a theory de- 
scribing the interaction of two stable particles A and B (or multiplets thereof) 
of different masses, itla and m#, respectively. One may find that the r.h.s. of 
Eq. (1) can be expressed by using only the second term if the integral path 
is extended to the complex q plane such as — oo + i( — » +oo + i( with 
m/2 < ( < m. In fact, the shift of the integral path back to the real go axis 
(— oo — ► +oo) picks up a residue associated with F{y) at the pole v = +vb, 
which leads to the first term of Eq. (1). In this sense, we may call the first 
term in Eq. (1) the pole term. Then, we may speculate that the asymptotic 
formula for the A-particle mass shift for the case vha ^> tub is written as 
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with Fab(v), the forward scattering amplitude for the A + B — ► A + B process, 
and in this case vb = m B /2mA- For instance, by identifying A with the nucleon 
and B with the pion, and exploiting the residue of Fn w (i>) at v — +ub, one 
would then reproduce the nucleon mass shift formula presented by Liischer in 
his Cargese lectures [9]. 

Among the applications of ChPT, the QCDSF-UKQCD collaboration esti- 
mated the finite volume effect on the nucleon mass from data in Nf = 2 
lattice QCD [6], applying the mass shift formula derived within ChPT. Along 
the way, they found however, when expressing their formula in terms of F^^{y) 
in the same order of ChPT, that the factor of the pole term is twice larger 
than that of Liischer's in Ref. [9], although the rest is the same. There seems 
to be no mistake in the formula in Ref. [6], at least within the infrared regu- 
larization scheme [1]. On the other hand, Liischer's formula being considered 
as general such that it can also be applicable to ChPT at any orders, this 
discrepancy poses a structural question. 



In this paper, we thus investigate the mass shift formula for the interacting two 
stable particle system along the lines of Liischer's proof for a self-interacting 
bosonic theory. We will not assume m^ ^> ttlb (as in the physical case above) 
because this would obscure the general structure of the formula, so that we 
will keep terms which would be relevant if the masses were close. In Sect. 2, we 
first describe the result on the generalized mass shift formulae for the boson 
and fermion in the boson-boson and fermion-boson systems, respectively. In 
Sect. 3, we then provide a part of the derivation by selecting a typical diagram 
for the self energy in the two particle system. We discuss the nucleon mass 
shift formula in the N-ir system in Sect. 4. A summary is given in Sect. 5. 
Notations used in this paper are summarized in Appendix A. 



2 Asymptotic formulae for finite size mass shift 



In this section, we describe the result on the generalized mass shift formulae 
for the boson and fermion in the boson-boson and fermion-boson systems, 
respectively. 

The physical mass of a stable particle is given by the position of the pole of the 
propagator of an asymptotic field. In the framework of perturbation theory the 
pole is shifted from the bare one due to the self energy arising from polarization 
effects from the virtual particles. In finite volume, the expressions for the 
self energy involve sums over discrete spatial loop momenta, q(L) = 27rn/L 
(n G Z 3 ), instead of integrals. By using the Poisson summation formula, such 
a summation can be rewritten again as an integral with another summation 
over integer vectors rh G Z 3 and an exponential factor depending on L: 
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where f(q) is a function composed of propagators and vertex functions. Then, 
the difference of the self energies between the finite and infinite volumes, ap- 
pearing in the definition of the mass shift as described below, can be defined 
by the sum over \m\ ^ 0, since \m\ — corresponds to the integral in infinite 
volume [11]. The asymptotic formula at large L, which is of our interest, is 
given by the contribution of \m\ = 1. (next leading contribution is \m\ = \/2.) 
In this case, the exponential factor in Eq. (4) is reduced to 2X)f=i cos(Lgj). 
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Fig. 1. Self-energy diagrams which contribute to the mass shift formula in the <f>A-<f>B 
system. A solid line with an empty circle corresponds to the propagator of 4>a and 
a dashed line to that of 4>b- A shaded blob is a vertex function, see Appendix A for 
the definition. It is assumed that 4>A carries a conserved charge. 

2.1 Boson mass shift formula in the 4>a-<Pb system 



As the finite and infinite volume masses Ma{L), itla of A-particle are defined 
from the poles of the propagators of 4>a (see Appendix A), Ami(L) = Ma(L) — 
ijia is reduced to 
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[E L (p) - E(p)] + 0((Am A ) 2 ) at p = (im A , 
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where ^l(p) and E(p) stand for the self energies of <f>A i n the finite and infinite 
volumes, respectively. We work in Euclidean space as in Ref. [10]. Here, since 
the fermion mass shift formula in the fermion-boson system is of our final 
interest, we assume that only <pA carries a conserved charge so that interaction 
induced only by the 3-point vertex AAB and 4-point charge conserving vertices 
are taken into account. In such a system, the difference of the self energy of 
4>a i n the finite and infinite volumes at any orders in perturbation theory is 
provided by the diagrams listed in Fig. 1. For tub jw, a £ (0, 1], Eq. (5) is 
reduced to 
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Fig. 2. Ingredients of F AA (y) in the (j) A -<p B system. The labels represent the corre- 
spondence with the self-energy diagrams in Fig. 1. The diagram (bl) in Fig. 1 is 
responsible for the two types of amplitudes. 
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where FabW) an d FaaW) denote the forward scattering amplitudes of the 
processes <f> A {p) + 4> B (q) -> <j> A (p) + <f> B (q) and <f) A (p) + <j> A (q) -»■ <f> A {p) + 
4>a{c() i n the infinite volume, respectively. The ingredients of these scattering 
amplitudes are graphically represented in Figs. 2 and 3. Xaab is an effective 
renormalized coupling denned from the residue of F AB (v) at v — ±u B = 
±m B /2m A as 
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The pole term in Eq. (6) is associated with the evaluation of the self-energy 
diagram (bl) in Fig. 1. The error term is defined by 



m> \jl(m\ - v B ) , 
which is due to the neglect of \m\ > y/2 contributions. If m B /m A < ol c 



2 — \/2 w 0.765, we may neglect the second term in Eq. (6), which decays 
more rapidly than the error term at large L. Note that the pole term in 
Eq. (6) is twice larger than that in Eq. (3) and two types of forward scattering 
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Fig. 3. Ingredients of Fab(v) in the ^A-^B system (or in the ^a-'/'b system with ar- 
rows to the solid lines). The labels represent the correspondence with the self-energy 
diagrams in Fig. 1. 

amplitudes contribute to the formula in Eq. (6). 

2.2 Fermion mass shift formula in the ^a-<Pb system 



From the poles of the propagators of ^a in the finite and infinite volumes (see 
Appendix A), we obtain 



Am A = u(p, r) [E L (p) - S(p)]w(p, r) + 0((Am A ) 2 ) at p = (im A , ) ,(9) 

2m A 

where u(p,r) and u(p,r) are spinors of ^ a with the spin r. The mass shift is 
then written down as 



Am A {L) 



X AAB c -L y /rn^-uj 

2u B 



871171 aL 

f 00 ^ e- L ^^^{F AA (iqo) + F AA (iq )} 

J-oo Z7T 



+ r d Ja e- L V^<FAB^ ) 

J-oo All 



0(e~ Lm ) , (10) 



where m is the same as in Eq. (8). Now, three types of the forward scattering 
amplitudes contribute to the formula, F AB {v): ty A (p) + 4>B(q) — > ^A(p) + (pB(q) 
(Comptontype), F AA (v)\ tf A {p) + tf A (q) -► ^ A (p) + ^A(q) (Miller type), and 
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Fig. 4. Self-energy diagrams which contribute to the mass shift formula in the ^ a-<Pb 
system. A solid arrow with an empty circle corresponds to the propagator of \I/ a- 
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Fig. 5. Ingredients of Faa(v) in the \P a-<Pb system. The labels represent the corre- 
spondence with the self-energy diagrams in Fig. 4. 



F AA (v): m A {p) + ^A{q) -+ y A (p) + ^A(q) (Bhabhatype) (see, Figs. 3, 5 and 6). 
Again we may neglect the second term if mB/rriA < ol c . Apart from the relative 
minus sign in front of the second term, which is due to Fermi statistics, the 
formula is almost the same as for the boson. The factor of the pole term is 
twice larger than that in Eq. (3). 
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Fig. 6. Ingredients of F AA (v) in the ^> a-^b system. The labels represent the corre- 
spondence with the self-energy diagrams in Fig. 4. 

3 Outline of derivation 



In this section, we provide a part of the derivation of the mass shift formulae 
in Eqs. (6) and (10) for the case 



a = ^ 6(0,1] 
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We shall here concentrate on the evaluation of the self-energy diagram (bl) in 
Figs. 1 and 4, which is typical in the two particle system and finally provides 
the pole term 0(e~ Lfl ) with 
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The definitions of the propagators and the vertex functions used in the follow- 
ing calculations are summarized in Appendix A. Other self-energy diagrams 
can be evaluated in a similar way. The results are collected in Appendix B. 



3. 1 Boson mass shift formula 



The self-energy diagram (bl) in Fig. 1 with \m\ — 1 (see Eq. (4)) is given by 
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p=(im A , 0) 



where 5 G (0, |] is a parameter chosen optimally below. We are working in 
Euclidean space so that % = (f = (<7o?<f ) and q = (qi,q2,Q3)- Inner product 
is defined by pq = poqo + p ■ q = irriAqo- Due to rotational invariance among 
q±, q 2 and q$ one can rewrite the cosine factor in Eq. (13) as 2 X^=i cos LQj ~^ 
6Ree iLqi . 

Firstly, we integrate out q±. To do this, we extend q\ to a complex variable and 
perform complex contour integration. Here, the integral path must be chosen 
appropriately by considering the analyticity properties of the integrand, the 
vertex functions and the propagators, in the complex q\ plane. 

A complex domain D = {(p, q) G C 4 x C 4 }, to which the vertex function 
Faab(TP, ±(1 — S)p ± q\ ±<5p T q) initially defined for (p,q) G M 4 x R 4 is 
analytically extended, can be found as follows. The basic assumption here is 
that the vertex function at any orders in perturbation theory consists of a set 
of A and B lines (free propagators). The /-th A or B line is then parametrized 
as [(k{l) + r(l)) 2 + mW' 1 or [(k(l) + r(l)) 2 + m|] _1 , where k(l) is the flow of 
external momentum given by a combination of complex variables p and q, and 
r(l) is a combination of internal loop momenta to be integrated out, which is 
a real variable in Euclidean space. Thus, we find that the vertex function has 
no singularity if (Im k(l)) 2 < m 2 A and (Im k(l)) 2 < m 2 B are satisfied for all A 
and B lines. In order to find the possible choices of k(l), we label the three 
bare vertices where the external momenta, ^fp, ±(1 — 5)p±q and ±Sp^fq, are 
plugged in as 01, a 2 and b, respectively (e.g. 01 = a 2 = b at the tree level). We 
consider here the case that A-particle carries a conserved charge. In this case, 
there always exists a set of A lines connecting a\ and a 2 - We can then choose 
k(l) up to overall sign as follows. If the connected A lines flow through b, such 
A lines carry (1 — 5)p — \q between a,\ and b, and (1 — 25)p + \q between b 
and a-i- The other lines then carry 5p + |g or 0. If the connected A lines do 
not flow through b, A lines carry (1 — S)p + |g between a,\ and a 2 , and the 
other lines carry Sp, |g, Sp — \q or 0. The analytic complex domain D is then 
specified by inserting these k{l) into the above two inequalities. Practically, it 
is enough to consider the cases 



Im{(l-5)p±^}) <m 2 A , (14) 

1 \ 2 
lm{5p±-q}\ <m|, (15) 

so that all other cases are fulfilled. 2 Solving these inequalities with Imp = 

2 In a similar way, one can discuss the analyticity properties of Taab for the case 
neither (f>A nor (f>B carry the conserved charge. In this case, the coupling ABB 
can be present and therefore the r.h.s. of inequality (14) must be replaced by m 2 B 
for a < 1. Then, 5 = 1/2 will be an optimal choice, which provides the analytic 



?7iA, P = 0, and (go, Q±) £ M 3 , where g_i_ = (52, 93), we find that the choice 



a 2 
5 = y (16) 

provides a maximum analytic strip of Im gi, in which the vertex function has 
no singularity: 



a 2 



< Imgi < 2m B dl- — = 2/x . (17) 

It is useful to find that if a = 1 (m^ = mg = m), the parametrization of 
external momenta (S = 1/2) and the analytic strip of Im q\ (maximum value 
is \/Zm) are reduced to the case as discussed in Ref. [10] for the identical 
bosonic system. 

Inserting Eq. (16) to the propagators Ga and Gb in Eq. (13), we find that 
these propagators possess poles of one-particle states at 
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in the complex q 1 upper half plane, respectively. 3 Note that these two poles 
happen to degenerate for q = as 



qr = q'i' = i^qi + ^- (22) 

To avoid the double pole in the complex q\ plane, we may deform the q§ 
integration path around q = into an infinitesimal half circle in the complex 
lower half plane (see Figs 7 and 8). For the purpose of controlling the error 
term associated with integration up to 0(e~ Lm ) (with in = V2fi), we consider 
go and q± in the ball 



strip < Im q\ < 2msyl — l/(4a 2 ). Clearly this strip can be defined only for 
a e (1/2,1]. 
3 For a, b G R, 



q = i\Ja + ib 



(Va 2 + b 2 - o)/2 + iJ(Va 2 + b 2 + o)/2 if b > , (20) 



+V (Va 2 + b 2 - o)/2 + ^(Va 2 + b 2 + o)/2 if 6 < . (21) 
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In this ball it turns out that Im q\ and Im q\ are inside the strips 
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Now, we choose a contour depending on a. li a G [a c , 1], where a c 



2 — y^ ps 0.765, we take a path which goes along the real q\ line and the line 
Im q 1 = J ji 2 + m 2 A closed at ±oo. In this case, both q[ and q[ are com- 



pletely inside the contour. If a G (0, a c ), we take the line Im q 1 = Jfi 2 + m 



'B 



instead of Jfj, 2 + m 2 A . In this case, while q[ is completely inside the contour, 
q[ locates partially outside the contour. However, since even the lower bound 
of Im q[ exceeds m for a < a c at next step of the complex contour integration 
of go, we may neglect the q[ contribution here. In both cases the contribution 
from the path where Im q\ ^ is negligible, since it is more rapidly decaying 
at large L than the error term 0(e~ Lm ). Then, from the residues at q[ and 
q\ , we obtain 
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where 



C=« 



jg» = M 



dqod 2 q ± e iLqi 
b (2vr) 3 2gi 

dq d 2 q± e iLqi 



Faab(~P, (1 - S)p + q;5p- q)G B (5p - q) 



xTaab(p, -(1 - S)p - q; -Sp + q)\ qi=q (A) 



Faab(~P, (1 - S)p + q;5p- q) 
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b (2tt) 3 2gi 
xG A ((l - 5)p + g)r^ B (p, -(1 - <J)p - g; -<fc + g)| m= „(B) . (2* 



As mentioned above, J fel is relevant only for the case a G [a c , 1]. Note that 
Gb(8p — q)\ _ (A) and GU((1 — <5)_p H- g) I _ <s) have a common pole at q = 0. 
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Fig. 7. qo integral contour for 1^ . 

Secondly, we perform the complex contour integration of go- Hi is evaluated 
along the contour in Fig. 7. The contributions from the paths BC and D'A 
are negligible at large L, since Im q[ > rh along these lines. The choice of the 
contour CD' (given by the shift go ~~ ► %~ *(1 — &) m A) is possible because both 

inequalities (14) and (15) are satisfied as long as a > y 3 — v7 ~ 0.595 and 
thus the vertex function has no singularity inside the contour for a G [a c , 1]. 
As there is no pole inside the contour, the integral can be written as that 
along the path CD'. The argument of the propagator Gb and q[ are then 
shifted as Sp 



q — > p — q and q\ 



(•4) 
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q\ + m A , respectively, where q 



satisfies the on-shell condition q = —vf\r A . The integral is then represented as 
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where 



f aa (iQo) = F A ab(-p, q;p- q)G B {p - g) 
xTaab(p, -g; -p + q) | P 2 =</2= _ 



(30) 



is a one-particle- irreducible (1PI) part of the forward scattering amplitude of 
Faa( v ) with v = pq/m A = iq (see the diagram (bl-2) in Fig. 2). 

( B) 

7 61 is then evaluated along the contour in Fig. 8. Again the contributions 

( R) 

from the paths BC and DA are negligible at large L, since Im q\ > rh along 
these lines. The integral along the path CD is parametrized by shifting the 
momentum as go — > go + iSrriA, which is possible because the vertex function 
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Fig. 8. go integral contour for i^ . 
is analytic inside the contour for the whole range of a G (0, 1]. The argument 



S B ) 



of G A and q\ are shifted as ( 1 — S)p + g — ► p + g and gj 



(s) 



_ gg + gi + m| 

and then g 2 = —m 2 B . Now, since the pole at go = is inside the contour, we 
also have a residue contribution. Thus the integral becomes 
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(32) 



is a 1PI part of the forward scattering amplitudes of F AB {y) (see the diagram 
(bl-1) in Fig. 3). Note that F AB '(v) has a pole at v — +u B . The coupling 
Xaab in the first term in Eq. (31) is the effective renormalized coupling defined 
through the relation 



X 2 

A AAB 



lim {v 2 -v B )F AB {y) = lim (v 2 - vl)F { tl' l) (u) 
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where all legs of the vertex function is on the mass shell, p 2 = —m\, {p + q) 2 = 
—m\ and q 2 = —m 2 B . The integral region of q± in the first term in Eq. (31) is 
defined by 



{q± e M 2 | q\ < /i 2 } . 



(34) 



One may find that the 1PI amplitudes corresponding to (bl-1) in Fig. 2 and 
(bl-2) in Fig. 3, which are necessary to build up the forward scattering am- 
plitudes in the mass shift formula, are apparently absent in 1^. However, 
we find that these amplitudes can be included by using the crossing rela- 
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F^ A ~ 2 \u) and F AB 
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ihi - 1) (-u) = F$- 2 \u); one can replace 

,(61-1), 



(61-2), 



FT A - Z \v) in Eq. (29) by (F^'(u)+F^'(u))/2 and F^\u) m Eq. (31) 

hj{Ft B - l \u) + Ft B - 2 \u))/2. 

Finally, we carry out the q± integration in Eqs. (29) and (31) using the formula 
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where the integral region is extended from B or B' to the infinity, because the 
additional contributions are always up to the order of the error term. Hence, 
we end up with 
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By combining the contributions from the other self-energy diagrams in Fig. 1 
(see Appendix B), we complete Eq. (6). 



3.2 Fermion mass shift formula 



Let us next consider the fermion mass shift formula in the ^a~4>b system. The 
self-energy diagram (bl) in Fig. 4 with \m\ = 1, sandwiched by the spinors 
u(p,r) and u(p,r), is given by 
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f d 4 q ( 3 \ _ 

Ibl = J J^kY 2 ^ COsL * u(p,r)T AAB (-p,(l-5)p + q;5p-q) 

xS A ((l - S)p + q)G B (5p - q)T AAB (p, -(1 - S)p - q; -5p + q) 



xu(p,r] 



(37) 

p=(iniA, 0) 



Since the pole structure of the integrand is the same as in the (f> A -(f) B system, 
we can carry out the momentum integration in a similar way. 

After several integration steps, we then arrive at the expressions (see Figs. 5 
and 3, respectively): 



#> = - 3 / y> e : LV ^%{F^) +iff faOHO(e-»») , 



(2^r) 3 2Jql + ql + m\ 



(38) 
^Ju^o) = -J2^(p,r)T a / AB (-p,q;p - q)llp(q, s)G B (p - q) 

s 

xu^q, s)T J J AB (p, -q; -p + q)us(p,r)\ p2=q 2 = _ m 2 A , (39) 

F AA( i( lo) = Yl ««(P> r ) r AAB(-P, -Q", P + Q)vp(q, s)G B {p + q) 

s 

xv^(q, s)T^ AAB (p,q; -p - q)u 5 (p,r)\ p 2 =q 2 = _ m 2 A , (40) 



and 



j(B) _ 3 / <*V e- L V^ X AAB 



'61 



(2tt) 2 2^ + ^ 2z/ B 



+ 3 / ^V '-T^ {f ffi-> (<gb ) + ^(fc), + O(e-) 

W 3 2 A 2 + gi + m| 



(41) 
F AB~ 1] ( i( lo) = u a {p,r)Vf AB (-p,p + q; -q)S^{p + g) 

>< r AAB(P,-p-(i;(l)MPi r )\p 2 =-m\,q2=-ml , (42) 

f ab~ 2) (iQo) = u a {p, r)T*f AB (-p,p - q; q)S ( A h (p - q) 

xT aab(p^ -P + T, -q)u s (p, r)\ p 2 = _ m 2 2= _ m 2 , (43) 



"A' * — "*B 



where we have used the relations m A — i</[ = X^ M (<?> s)u(q, s) and m^ + i^ = 
— S s ' i; (Q') s )v(q, s) in Eqs. (39) and (40), respectively. Repeated Greek letters, 
corresponding to spinor components, are implicitly summed. The overall minus 
sign in F AA is due to Fermi statistics. The effective coupling \ AA b is defined 
as in Eq. (33). The advantage of such a definition of the coupling in the 
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fermion-boson system is that the model dependent matrix structure of the 
vertex function, such as {/, 7^, a^, 757^,75}, is involved in Xaab through the 
definition of Fab- 

After integrating out q±, we obtain the expression 



'61 



4ttL 



X 2 

A AAB -L\l 

2u B 



^ e -L^/gf 

, 2tt 



; {F { A b l\iqo) + Ffl\iqo)} 



J -00 2% 



(61-1) 



(*9o 



? (61-2) 
AS 



(*?o)} 



0(e 



—Lm\ 



.(44) 



We find that the expression is almost equivalent to Eq. (36) apart from the 
relative minus sign in the second term. By combining the contribution from 
the other self-energy diagrams (see Appendix B), we complete Eq. (10). 



4 The nucleon mass shift 



As an application of the fermion mass shift formula in Eq. (10), let us discuss 
the nucleon mass shift in the N-n system. Since the formula is expected to 
hold nonperturbatively, it is interesting to estimate the mass shift by inserting 
the N-7i scattering amplitude which is known from experiment. Therefore, the 
following analysis can be regarded as an estimate of the finite volume effect 
on the nucleon mass when the realistic pion mass is achieved in lattice QCD 
simulations with dynamical fermions. 

According to Holder [12], the subthreshold expansion of the N-tt forward scat- 
tering amplitude around v = is parametrized as 



F N n(v) = Qm N D + (u) 
where 



(45) 



D + (u) 



9' 



vl 



m N v% 



+ d+ 



00 



m„ 



df m;V + d+ m;V + 0(u 6 ) . (46) 



The isospin sum is taken in Eq. (45) and the effect of isospin symmetry 
breaking is neglected. The masses of the nucleon and pion are then given by 
rriN = 938 MeV and m^ = 140 MeV. The coupling constant is g 2 /47r = 14.3. 
The first term in Eq. (46) is identified with the pseudovector nucleon Born 
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term with vb = m^/2mN ~ 0.07771^. The coefficients of the other terms are 
given by d& = -1.46(10), dfo = 1.12(2) and dt = 0.200(5) [12]. We only take 
into account the mean of these values hereafter. Note that although Eq. (45) 
is defined in Minkowski space-time such that v is a real variable here, it can 
directly be inserted to Eq. (10) by rotating this variable as v = iq®, where go is 
the integral variable in Eq. (10). The effective coupling is computed by using 
Eq. (7) as 



Mm (z, 2 - ul)F Nn (u) = 3mjL = _ 6 ^| . (47) 

v—*±ub Z 

Since m n /mN = 0.149 < a c , we may neglect the second term in Eq. (10). It is 
clear that since the factor of the pole term in Eq. (10) is already twice larger 
than that in Eq. (3), we obtain a different expression from that in Ref. [9]. In 
fact, the mass shift formula, divided by the nucleon mass itself, is reduced to 



5(£, = Lm n ) = AmN(L)/mN 




m N J 

2 



^-) r dy e-t^F N «(im«y) + 0{e 



Lm\ 



167r 2 £ \rriN ■ 
= MO + MO + MO + 0{e~ Lm ) , (48) 

where <M0 ls the P°l e t erm 5 an d <M0 an d $r(Q correspond to the contri- 
butions of the pseudovector Born term and the rest in Eq. (46). Clearly, this 
formula recovers the expression given in the framework of ChPT [6]. 

We plot <5(£) in Fig. 9, where £ = 1 corresponds to L — 1.4 fm. We find that 
5(C,) suffers strongly from higher order contributions of v in the range £ < 1. 
For instance, <M0 causes the negative mass shift within the leading mass 
shift formula (\m\ = 1). In this range, the contribution from \m\ > a/2 to the 
formula, of course, will not be negligible. On the other hand, 5(£) seems to be 
mostly described by <M0 as £ increases. However we notice that this is due 
to the cancellation between (MO an d <M0- 

Finally, let us discuss the relation between the general fermion mass shift 
formula in Eq. (10) and the nucleon mass shift formula derived in Ref. [6]. 
The diagrams that they evaluated within ChPT correspond to (bl) and 
(a2) in Fig. 4, where the vertex functions were treated as certain coupling 
constants. The diagram (bl) was evaluated by introducing the Feynman 
parameter x, where its initial integral range x G [0, 1] was modified as 
[0, 1] = (— oo, oo) — (— oo, 0) — (1, oo) and the terms of 0(e~ LniN ) which mainly 
originate from the integral in the interval (l,oo) were omitted according to 
the infrared regularization [1]. Then, from the integrals in the range (— oo, oo) 
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Fig. 9. The nucleoli mass shift as a function of £ = Lm^. 

and (—oo,0), they obtained the terms corresponding to <5p(£) and <5b(£) in 
the mass shift formula, respectively. The evaluation of the diagram (a2) was 
straightforward and they obtained terms for 8r(£). In our point of view this is 
reasonable since the diagram (a2) has nothing to do with the pole term. Here, 
it is interesting to notice that the omitted terms according to the infrared 
regularization can exactly be used to reconstruct the terms involving Fnn(v) 
and F N ft(v), which we have also omitted in Eq. (48) because of a < a c . In 
this sense, the general formula in Eq. (10) covers all possible contributions to 
the mass shift, some of which may be neglected depending on the physical 
situation. 



5 Summary 



We have studied the general finite size mass shift formula for the two stable 
distinguishable particle system in a periodic finite volume along the lines of 
Luscher's proof for an identical bosonic theory. The main results are Eqs. (6) 
and (10). What we have clarified are the following. 

The boson mass shift formula for <J)a in the 4>a~4>b system is related to two types 
of the forward scattering amplitudes, Fab(v) and Faa(v). If m B /mA < a c = 

y 2 — \f2 m 0.765, the contribution from the Faa{v) term becomes smaller 
than the order of the error term. It has turned out that the speculation such 
as in Eq. (3) is not valid for a two particle system. However, the mechanism 
how the mass shift, proportional to the difference of the self energy S^ — E, 



is related to the forward scattering amplitude is the same as in the identical 
bosonic system [10]. In principle, the difference of the self energy consists of the 
propagators of 0^ and 0#, one of which is broken into external legs on the mass 
shell after integration over the loop momentum. Then, a self-energy diagram is 
related to one (or some) of the 1PI parts of the forward scattering amplitude. 
It is to be noted that in the two particle system, unless the crossing symmetry 
among the 1PI parts of the amplitude is made manifest, some of the 1PI 
parts necessary to build up the forward scattering amplitudes will apparently 
be missed although such an expression is mathematically consistent. (In the 
identical particle system, the crossing symmetry is automatically manifest.) 

The fermion mass shift formula for ^a in the \1/ a-<Pb system is related to three 
types of the forward scattering amplitudes, Fab(v), Faa(v) and F AA -{y). If 
'mB/mA < « c , the contributions from the Faa(v) and F AA {y) terms become 
smaller than the order of the error term. Again, crossing symmetry must be 
taken into account to relate the self-energy diagrams to all 1PI parts of the 
forward scattering amplitudes. When the fermion propagator is broken into 
external legs on the mass shell, the relative minus sign appears in front of the 
corresponding terms. 

In both cases, Xaab = unless the self energy contains the (bl) diagram, 
since it means that Fab{v) has no pole. 

Finally, we have written down the nucleon mass shift formula by applying the 
fermion mass shift formula in Eq. (10) to the N-tt system, thereby we have 
found that the pole term is underestimated by factor two in Liischer's formula 
in Ref. [9] and the formula derived in ChPT [6] is reproduced. 

There are now general finite size mass shift formulae for the two particle 
systems which are valid to all orders in perturbation theory within \m\ = 1, in 
addition to that in the identical bosonic system [10]. For every case all these 
formulae are obtained in the same way by carefully analyzing the appropriate 
set of self-energy diagrams. 
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A Notation 

In this appendix, we summarize our notations. 

A.l Boson propagator 

The boson propagator in infinite volume in Euclidean space is defined by 



G{x-y) = J-^G{p)e^-y\ (A.l) 

°(P) = 2 ^ I w , » (A.2) 

p A + m z — h{p) 

E(p) = ^E(p)=0 for p 2 = -m\ (A.3) 

where m is the physical mass. The renormalization conditions (A.3) mean that 
the propagator has pole at p 2 = —m 2 with unit residue. We assume that there 
are no bound states so that this pole is the only singularity of G(p) below the 
two-particle threshold. 

In the finite volume of the size L 3 , by imposing periodic boundary conditions 
for the spatial directions, we have 



G L {x-y) = ^Y.j d ^G L W ip(x - y) , (A.4) 

p 

G L(P) = 2^_ 2 V f \ ' ( A - 5 ) 

p2 _|_ m Z _ Li(pJ 

where the momenta are quantized as p = 2irn/L (n G Z 3 ). The summation 
of discrete momenta can be reformulated as the integral by using the Poisson 
summation formula (see Eq. (4)). We expect that the pole position of Gl{p) 
is shifted from p 2 = —m 2 to p 2 = —M(L) 2 with M(L) = m + Am(L). 

In the 4>a-4>b system, by expanding Sl(p) in powers of Am,A{L), we then solve 



G L (p) 1 \ p 2 = _ Ma{l) 2 = , (A. 6) 

and obtain Eq. (5) for Ai71a{L). 
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A. 2 Fermion propagator 

The fermion propagator in infinite volume in Euclidean space is 



S(x -y)=J -0^S(p)e^^ , (A.7) 

S(p) = ... l w s , (A.8) 

d 

^(p) = tt7^(p) = f° r ip = ~ m > (A-9) 

op 

where p" = JqPo + 7j S £>i- Euclidean Dirac matrices are defined by 7jf = 70 and 
yf = i'ji = —iY (i = 1,2,3), which satisfy {7^,7^} = 25^ v . As similar to 
the boson propagator, we also assume that there are no bound states so that 
the pole at if 1 = —m is the only singularity of S(p) below the two-particle 
threshold. 

In the finite volume of the size L 3 with periodic boundary conditions, we have 



S L (x-y) = y^E J '' ^S L {p)e^ x - y) , (A.10) 

v 

Sl ^ = -wt-L. X v n ' (A - n) 

if + m — I-* l{P) 

where we expect that the pole position is shifted from if 1 = —m to if 1 = —M(L) 
with M(L) =m + Am(L). 

In the ty a-'Pb system, we then solve the equation 



u(p,r)S L (p) l u(p,r)\ p 2 = _ MA(L) 2 =0 , (A.12) 

where u(p,r) and u(p,r) are spinors of ^a defined as the solution of the free 
Dirac equation 



{if 1 + m)u(p, r) = , u(p, r) (ip" + m) — , 

u(p, r)u(p, s) = 2m5 rs , (A. 13) 

and obtain Eq. (9) for Ai71a{L). 
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Fig. A.l. Vertex functions used in Figs. 1~6. 
v4 . 5 Vertex functions 



In Fig. A.l, we list all types of the vertex functions used in Figs. 1~6. Arrows 
represent the direction of momentum flow. Outgoing momenta from the vertex 
are parametrized to have positive sign. Greek letters denote spinor indices for 
the ^ a~4>b system. 



B Ingredients of the formulae 



We list all contributions to £z,(p) — S(p) in the (pA-rf'B system (Fig. 1) and 
u(p, r)[Ei(p) — T,(p)]u(p, r) in the \I/ a-4>b system (Fig. 4). 



B.l 4>a-4>b system 



'ol 



dto-Ly/. 



{al). 



e-^^m lF ^ >{iqo)+0{e 

47T.L J-oo Z7T 

F AA(iQo) = r AAAA{-p,P,-q,q)\p2= q 2=-ml 



—Lm\ 



(B.l) 



l a2 



A-kL 



dto-Ly/. 



2tt 



$+»% Fffliith) + Ok-™) 



(B.2) 
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^ ' / 1 /i 1 II ,.,,,; t I ) I • _ . / / t 1 II n 9 9 9 

n z — — Tn^ n-^ — — -m-^ 



^Ab (m) = r A ABB(-p,p;-q,q)\ 



p*=—m., q^=-m^ B ) 



<bl 



AuL 



X AAB -Ljml-uj 
1VB 



+ f ^-Vra {f fi-i) (jft ) + jjj-^j 



-OO 

/■OO 



f aa~ 1] '(*9o) = r^Bt-p, -g;p + g)G B (p + g) 

xr j4j4iJ (p,g;-p-g)| p 2 = g2 = _ m 2 i , 

Fab'^ (ilo) = F A ab{-P,P + q; ~q)G A {p + q) 

xTaab(p,-p- q;q)\p2=-m 2 A ,q2=-m% , 

FjtiT^iiqo) = Faab{-P,P- q; q)G A {p - q) 

\p 2 =—m., q- L =—m,g i 



+ 0(e~ Lm ), (B.3) 



xTaab(p, —p + q\ — g)lp 2 ^-»<- „j_ ._,„- 



/— OO 

F i c ]W) = r^ B (-p,p;0)G ! B (o)r j4v4s (-g,g;0)| P 2 =9 2 = _ m 2 , 



/cl = 4^L loo f^ V?0+mi nf(^o) + 0(e-^) , (B.4) 



'c2 



r ^e^v^^l Fi^(igo) + 0(e~ L n , (B.5) 

J-oo ZlT 



AtiL 

F ab{Wo) = T A AB{-p,p;0)G B (0)T BBB (0,-q,q)\ r ._ _,„-,,, _,„- 



The forward scattering amplitudes are defined by 

i^(»«b) = F^iiqo) + F^iiqo) + F A b \~ 2) \iq ) + F$(iq ) , (B.6) 
F AB (iqo) = F { A f(iq ) + F^ (iq ) + F^fo) + F$(iq ) . (B.7) 

B.2 ^a-'Pb system 



4i = ~ /_~ ^e^v^^jFif^o) + ^(^0)}+ 0(e"^) , (B.8) 

,(ol) 



F AA (*9o) = -^MP> r )MP> r ) M 7(<?, s)w 5 (<?, S) 
s 

xI "aaaa(-PiP, -g,g)|«2=«2=- m 2. 



Ip^=(jr' a =— m^ 5 
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2 =— m 2 A i 



FaA^Io) = X^a(p,r> / 3(P> r K(?> s)u s (q, s) 

s 

^ a AAAA(-v,VAi -g)| P 2= g 2=- 

4 2 = ^T r ^e- L V^F< F^\iq ) + 0(e- L n , (B.9) 

4:71 L J-oo Z7T 

Fab = Ua{p,r)T a / ABB {-p,p; -q,q)u p {p,r)\ p2= _ m 2 A!q2= ^ m 2 B , 



j _ 3 



X AAB Ljm'j-ul 

Ivb 



/oc 
-oo 



^e-Vi^ { irCS)( ift))+i r(«J) (jft))} 
+ r ^e-v/^l {FiT^o) + *- 2) (^o)}l + 0(e—) ,(B.10) 
F AA(iqo) = -J2u a (p,r)rf AB (-p,-q;p + q)up(q,s)G B (p + q) 

s 

XM 7 (g, s)r^4 B (p, g; -p - q)u s (p, r)\ p 2 =q 2 = _ m 2 A , 

F AA(iQo) = ^u a (p,r)rf AB (-p, q;p- q) v f3(q, s)G B (p - q) 

s 

xv^q, s)r J J AB (p, -q; -p + q)u$(p,r)\ p 2 =q 2 = _ m 2 A , 
F AB~ 1] ( i( lo) = u a (p,r)rf AB (-p,p + q; -q)S P A '{p + q) 

xT AAB(P^-p-Q^) u s(P^ r )\p 2 =^m% q 2 = - m % , 

F AB~ 2) (m) = u a (p,r)rf AB (-p,p - q; q)S A h (p - q) 

X T AAB (P, ~P + <?; ~q)us (P, r) \p>=- m \, q 2=-m% , 

/ci = ~ /_" |VV^i { FJ3> (<«,) + FgW)} + 0(e"^) , (B.ll) 

s 

xG B {^)u 1 {q,s)T 1 AAB (-q,q]Q)u & (p,r)\ p 2 =q 2 = _ m 2 A , 

F aaW) = M Q (p,^)r^ B (-p,p; 0)vp(q, s) 

xG B (0)v 7 (q, s)T A 6 AB (-q, q;0)u d (p,r)\ p 2 =q 2 = _ m 2 A , 

I c2 = -^ T r ^v^lfgW+Ote^) , (B.12) 

F AB{m) = u a (p,r)rf AB (-p,p;0)up(p,r)G B (0) 

xT B Bb(0, ~q, q)\ p 2=-m 2 A , q 2 =-m% ■ 

The forward scattering amplitudes are defined by 
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F AA (iq ) = F^\iq ) + F^\iq ) + FJ$(iq ) , (B.13) 

F AA (iqo) = F^\iq ) + F^\iq ) + F^(iq ) , (B.14) 

F AB (iqo) = F A f(zq ) + F A bl B - l \iq ) + F^~ 2) (iq ) + F^{iq Q ) . (B.15) 
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